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Abstract. 

We present the results of work in which we derived the 0(p 2 ) corrections to the Friedmann 
equations in the Randall-Sundrum I model. The effects of Golberger-Wise stabilization are taken 
into account. We surprisingly find that in the cases of inflation and radiation domination, the leading 
corrections on a given brane come exclusively from the effects of energy density located on the 
opposite brane. 



INTRODUCTION 

In recent years, there has been considerable interest in the idea that our world might be 
a 3-brane embedded in a higher dimensional bulk [1, 2, 3]. The fact that gravity alone is 
allowed to propagate through the bulk not only accounts for why the extra dimension(s) 
have so far avoided detection, but also provides an attractive solution to the infamous 
hierarchy problem. 

We will be concerned here mainly with the Randall-Sundrum I model, where two 
3-branes bound a slice of AdS(5) space. In this model, the weakness of gravity on our 
brane, the so-called "TeV brane" which has negative tension, is a consequence of the 
fact that the extra dimension is warped, i.e., the metric along different slices of the bulk 
parallel to the 3-branes will depend on their position along the extra dimension. 

An interesting feature of models with 3-branes in a 5-dimensional bulk is the fact 
that the Friedmann equations contain terms of higher order in the energy density than in 
standard cosmology [4, 5, 6, 7, 8]. More specifically, H 2 = + 0(p 2 ). These new 
terms could have important implications during inflation [9, 10, 11] and electroweak 
baryogenesis [12, 13]. Until now however, their exact form had only been worked out 
in the Randall-Sundrum II model, which doesn't solve the hierarchy problem. (In RS II, 
the extra dimension is infinitely large, and there is only a single 3-brane present). 

We present here the results of work [14] (see also [15]) whose aim was to find the 
0(p 2 ) corrections to the Hubble rate in the RS I model, taking into account the effect of 
the Goldberger-Wise (GW) mechanism [16] for stabilizing the extra dimension. 



EINSTEIN EQUATIONS 



The Einstein equations follow from the action 

S = Jd 5 x^^-^R-A+±d m <t>d m <t>-V(<i>)^ 

+ J d 4 x^g(j? mfi -V (<i>)) \ y=Q + J \y=V CD 

where K 2 = M~ 3 , and the potential V(4>) is left unspecified for now. The brane contri- 
butions are a sum of matter, represented by 

^,o~p*; ^ m ,i-p (2) 

and tension, which is the value of the brane's scalar potential V-(4>(y ; -)). The matter 
Lagrangians cannot be written explicitly for cosmological fluids, but their effect on the 
Einstein equations is specified through their stress-energy tensors (5). Our ansatz for the 
5-D metric has the form 

ds 2 = n 2 (t,y)dt 2 - a 2 (t ,y)Y^dx 2 -b 2 (t,y)dy 2 

i 

= e- 2N ^dt 2 -a (t) 2 e- 2A ^^dx 2 -b 2 (t,y)dy 2 . (3) 

i 

We will make a perturbative expansion in the energy densities p , p* of the branes around 
the static solution, where p = p* = 0: 

N(t,y) = A (v) + 5iV 1 (?,v) + 5^ 2 (?,v); A(t ,y) = A Q (y) + 8A^t ,y) + 8A 2 (t ,y) 
H*,y) = b + 8b l (t 1 y) + Sb 2 (t,y); <t>(t,y) =& (y) + 8& l (t,y) + 8& 2 (t,y). 

(4) 

The subscripts on the perturbations indicate their order in powers of p or p*. This 
ansatz is to be substituted into the scalar field equation and into the Einstein equations, 
G mn = K 2 T mn . (Since the scalar field equation can be derived from a combination of the 
Einstein equations, we will not worry about it any further). The stress energy tensor is 
T mn = gmn(V(<t>) + A) + d m ®d n ® - \d l ^d^g mn in the bulk. On the branes, T m n is given 
by 

V = ^Ldi a g(V + p,,V -p,,V -p,,V -p,,0) 

+ ^ T ^diag(y 1 +p,y 1 - jP ,y 1 -p,y 1 -p,o) (5) 

(Later we will assume that the potentials V Q and V l are very stiff and are vanishing at 
their minima, so they can be neglected.) 



It will be useful to define the following variables, which appear naturally in the Israel 
junction conditions: 



* 2 = 8A 2 -A^ - ^ 5* 2 - £ (s^ + ^ ^j T 2 = 8N' 2 - 8A 2 . 
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The analogous quantities at first order were found to be 
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in [15], and they satisfy the same boundary conditions as in (7,8), but with the replace- 
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In terms of these variables, we can write the second order Einstein equations as 
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where all the dependence on first order quantities squared is contained in the functions 
<^vy, <&y, ^05 ant * ^55- H ere ' me equivalent first order equations can be obtained by 
simply replacing the subscripts 2 for 1, and leaving out the functions 

We thus have at each order in the perturbative expansion a set of first order differential 
equations which, when combined with the boundary conditions, allows us to solve 
for the unknown functions T n and 8b n /b Q . (We will work in a gauge where the 
fluctuations 8<& n vanish). 

One final note regarding the physical value of the Hubble rate. Since the (00) com- 
ponent of the metric receives corrections when carrying out our perturbative expansion, 
we should use the time coordinate 



dx = 



n(t,l) 
n (t,l) 



dt = e SN 1 (t,l)-8N 2 (tA)-- dt 



(14) 



to define the Hubble rate. This means that physical Hubble rate is given by 

H^il + dN^ (15) 

rather than simply H = | as we might naively have expected. (See [14] for a more 
thourough discussion of this issue.) 



FRIEDMANN EQUATIONS 

Taking into account all that has been said in the previous section, we are now able to 
give the Friedmann equations to second order in p and p*: 



H 
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\y=l ~ — 
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dH 
~d% 



\y=0 



+4(l-3w,)(4 + 3w,)a 2 p* 2 + 4(l-3w)(4 + 3w)pp,))^ (16) 
P + P* + 3^2 ( 9 ( ! -3w,)(l + w,)p, 2 a 4 -(l-3(0)(7 + 3(0)p 2 
+2a 2 pp,[2(l-3(0)(2 + 3w)-3(l-3(0,)(l + w)] (17) 



dH 
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= -47rG( p(l + fi>)+p*(l + fife) 

y=l V 

+ (^ 2 (1 + -3w,)(13 + 9w,)p 2 + 9(l-3w)(l + w) 2 p 2 

(l-3w)(2(l + (0,)+2(l + (o)+6(l + (o) 2 + 3(l + co)(l + w,))pp,) ^ 

(18) 

= -47z;G(p(l + <0)+p*(l + <0*) 

+ T^S? f-4(l + ©)(!- 3(0)p 2 + 9(1 + co*) 2 (l - 3w*)C2 4 p* 



3m 2 Q 2 

+ a 2 [6(l-3©)(l + «) 2 + 2(l-3©)(l + ©)-2(l-3«)(l + <»*) 

4(1 -3cd*)(1 + ©) + 3(1 + ©*)(1 - 3fi>)(l + fi>) ]p*p ^ (19) 

where we have used 



p = -3H(p+p) = -3H(l + 0))p; SnG = k 2 (lb Q J^e'^dy^j ; 



(20) 



-3#(p* + p*) = -3#(l + (0*)p*; m 2 ^ ^fc 2 v^e 2 fc 2 a 2+2e ; (21) 



a = e -My =1 ). 



(22) 



It should also be noted that in order to obtain analytical results, we had to expand our 
expressions in powers of Q., and we have only kept the dominant terms. 

The careful reader will notice that the second order corrections vanish identically for 
radiation-like equations of state (&) = &)* = 1/3). In that case, we need to look at next to 
leading order terms to find non- vanishing corrections: 

A-. 

H \=o 

and the equations for dH/dx have vanishing corrections at this order in p and p*. 

DISCUSSION 

Let's now look at what our results mean. First, we note that the brane on which the 
hierarchy problem is solved is the one located at y = 1. It is thus natural to assume this 
is the brane we are living on. If we look strictly at the terms linear in p,p*, we see that 
in order for the current Hubble rate on our brane not to be completely dominated by the 
energy density on the other brane, we must assume that p* is presently very small or 
vanishing. 

If we now look at the second order terms on our brane during inflation (CO = — 1), 
we see that the only corrections come from terms involving p*. So unless we can come 
up with a model where p* was large in the past and becomes negligible in the present 
epoch, there can have been no appreciable effect coming from the second order terms 
during inflation. 

The situation during radiation domination is similar. The leading second order correc- 
tions on each brane involve the other brane's energy density. Furthermore, even if for the 
sake of argument we allow p* to have been present during this era, we know that it can't 
be more than 10% of the value of p by the time of nucleosynthesis. This means that 
we can ignore the Q^p* term, and that the term that is left over suppresses the Hubble 
rate. But in order to make the sphalerons go out of equilibrium at the electroweak phase 
transition, we would have needed to make the Hubble rate larger [12]. So once again, we 
find that the second order corrections do not seem to have any constructive applications 
in this model. 



CONCLUSION 

We have derived the second order corrections to the Friedmann equations in the Randall- 
Sundrum I model, taking into account the effect of the GW stabilization mechanism. We 
have found that the corrections during radiation domination have the wrong sign to be 
useful for electroweak baryogenesis. We have also found that any effect on a brane 



= —( (P* + P)+ P*-P)P* J (23) 

SnG ( , , 2kG ,~ , \ 
= — \b* + P)-tf^&Q*-P)p) (24) 



with an inflationary equation of state can only be coming from the opposite brane. Our 
approach being perturbative in nature, we must however keep in mind that our results 
should not be expected to hold when the energy density becomes greater than some 
critical density, which on the TeV brane is the (TeV) 4 scale. 



Acknowledgements 

We wish to thank Hassan Firouzjahi for enlightening discussions. J.V. is supported in 
part by a grant from Canada's NSERC. 



REFERENCES 

1. Arkani-Hamed, N.,Dimopoulos, S., and Dvali, G., Phys. Lett. B 429, 263 (1998) [hep-ph/9803315]. 
Antoniadis, I., Arkani-Hamed, N., Dimopoulos, S. and Dvali, G., Phys. Lett. B 436, 257 (1998) 
[hep-ph/9804398]. 

Arkani-Hamed, N.,Dimopoulos, S., and Dvali, G, Phys. Rev. D 59, 086004 (1999) [hep- 
ph/9807344]. 

2. Randall, L. and Sundrum, R., Phys. Rev. Lett. 83, 3370 (1999) [hep-ph/9905221]. 

3. Randall, L. and Sundrum, R., Phys. Rev. Lett. 83, 4690 (1999) [hep-th/9906064] . 

4. Binetruy, P., Deffayet, C. and Langlois, D., Nucl. Phys. B 565, 269 (2000) [hep-th/9905012]. 

5. Csaki, C, Graesser, M., Kolda, C. and Terning, J., Phys. Lett. B 462, 34 (1999) [hep-ph/9906513]. 

6. Cline, J. M., Grojean, C. and Servant, G, Phys. Rev. Lett. 83, 4245 (1999) [hep-ph/9906523]. 

7. Kanti, P., Kogan, I. I., Olive, K. A. and Pospelov, M., Phys. Lett. B 468, 31 (1999) [hep- 
ph/9909481]. 

8. Csaki, C, Graesser, M., Randall, L. and Terning, J., Phys. Rev. D 62, 045015 (2000) [hep- 
ph/99 11406]. 

9. Kaloper, N. and Linde, A. D., Phys. Rev. D 59, 101303 (1999) [hep-th/981 1 141]; 
Arkani-Hamed, N., Dimopoulos, S., Kaloper, N. and March-Russell, J., Nucl. Phys. B 567, 189 
(2000) [hep-ph/9903224]. Cline, J. M., Phys. Rev. D 61, 023513 (2000) [hep-ph/9904495]; 
Nihei, T., Phys. Lett. B 465, 81 (1999) [hep-ph/9905487]; 

Himemoto, Y. and Sasaki, M., Phys. Rev. D 63, 044015 (2001) [gr-qc/0010035]. Lukas, A. and 
Skinner, D., hep-th/0106190. 

10. Maartens, R., Wands, D., Bassett, B. A. and Heard, I., Phys. Rev. D 62, 041301 (2000) [hep- 
ph/99 12464]. 

11. Copeland, E. J., Liddle, A. R. and Lidsey, J. E., Phys. Rev. D 64, 023509 (2001) [astro-ph/0006421]. 

12. Servant, G, [hep-ph/0 11 2209]. 

13. Csaki, C, Graesser, M. L. and Kribs, G. D., Phys. Rev. D 63, 065002 (2001) [hep-th/0008151]. 

14. Cline, J.M., and Vinet, J. JHEP 02 (2002) 042 [hep-th/0201041]. 

15. Cline, J. M. and Firouzjahi, H., Phys. Lett. B 495, 271 (2000) [hep-th/0008185]. 

16. Goldberger, W. D. and Wise, M. B., Phys. Rev. Lett. 83, 4922 (1999) [hep-ph/9907447]. 



